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As the era of high precision cosmology approaches, the empirically determined power spectrnm 
of the microwave background anisotropy, Ci, will provide a crucial test for cosmological theories. 

We present a unified semi-analytic framework for the study of the statistical properties of the Ci 
coefficients computed from the results of balloon, ground based, and satellite experiments. An 
illustrative application shows that commonly used approximations bias the estimation of the baryon 
parameter Qb at the 1% level even for a satellite capturing as much as ~ 70% of the sky. 


During the next few wars ground based observations 
and balloon missions Q as well as satellite observa¬ 
tions promise exquisite determinations of the power- 
spectrum, Cl, of cosmic microwave baclground (CMB) 
anisotropies. Recently, pioneering work |g|-y has shown, 
in the context of inflationary cosmogonies, the tremen¬ 
dous impact these missions will have on our knowledge 
of cosmological parameters. It is therefore of paramount 
importance to study the statistical properties of these 
quantities in detail. The fact that the full solution of 
the associated likelihood problem poses severe compu¬ 
tational difficulties, has prompted many practical anal¬ 
yses to invoke one or more of the following simplifying 
assumptions: 1 ) the observed Ci are approximately inde¬ 
pendent due to nearly full and uniform sky coverage, 2 ) 
their sampling distributions do not change appreciably 
from the distributions which apply for the full sky, 
apart from a rescaling of the mean by the sky fraction 
to account for lost power, and 3) the sampling distri¬ 
butions are well-approximated by a Gaussian which has 
the same first and second moments as these rescaled 
distributions. 

There are several reasons to assess these approxima¬ 
tions and, if necessary, go beyond them. One reason is 
that in the short term balloon and ground based experi¬ 
ments will provide the leading edge science results to the 
field. Due to practical limitations we cannot hope to even 
come close to full sky coverage with these types of exper¬ 
iments. The observation regions are often ring-shaped or 
cover a circular region of the sky which subtends a small 
solid angle. It is an urgent matter to assess statistically 
how these imminent experiments can constrain the power 
spectrum and cosmological parameters. 

Further, the planned satellite missions are aiming to 
determine the Ci to sub-percentage accuracy. In the case 
of the Planck Surveyor mission, this has given us the hope 
of detecting small effects such as secondary anisotropies 
which are due to nonlinear gravitational effects on the 
CMB photons during the free-streaming epoch. This is 
an important issue because a detection would break oth¬ 
erwise present parameter degeneracies and allow a con¬ 


sistent parameter estimation from CMB data alone | 8 j. 

More generally, the impact CMB observations will have 
on cosmology makes it important to use approximations 
in a controlled way. For example, in the analysis of 
COBE-DMR data it was realized that using Gaussian 
approximations for quadratic quantities introduced sys¬ 
tematic biases [H. At the same time, we need approxi¬ 
mations to make feasible the analysis of the huge CMB 
data sets we expect in the coming years. 

To provide a quantitative basis for this discussion we 
present in this Letter a semi-analytic framework for the 
calculation of the sample statistics of the C/ for theories 
which generate a Gaussian CMB sky. This framework 
is exact for survey geometries and noise patterns which 
obey rotational symmetry about one, arbitrary axis, such 
as polar cap shaped regions, Galactic cuts, rings, and 
annuli of any size and at any latitude. We find that our 
methods also allow dealing with arbitrary noise patterns 
to very good accuracy and conclude by illustrating their 
use in a first application. 

The full sky of CMB temperature fluctuations can be 
expanded in spherical harmonics, Y/m, as|^ 

Til) = X X ^ImClmil) (1) 

i=0 rn 

where 7 denotes a unit vector pointing at polar angle 9 
and azimuth (j>. A Gaussian cosmological theory states 
that the aim are Gaussian distributed with zero mean 
and specified variance = {\aim\'^)- Hence, for 

noiseless, full sky measurements, each measured Ci inde¬ 
pendently follows a x^-distribution with 2 / -|- 1 degrees 
of freedom and mean 

Owing to Galactic foregrounds, limited surveying time 
or other constraints inherent in the experimental setup, 


*We assume that there is insignificant signal power in modes 
with I > Imax and use the convention that sums over m run 
from —Imax to Imax and all quantities with index Im vanish 
for m > 1. 
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the temperature map that comes out of an actual mea¬ 
surement will be incomplete. In addition, a given scan¬ 
ning strategy will produce a noise template. We model 
the noise as a Gaussian field T^v with zero mean which 
is independent from pixel to pixel and modulated by a 
spatially varying rms amplitude Wn{pi). Therefore the 
observed temperature anisotropy map is in fact 

fin) = Wij) [T( 7 ) + WNh)TNh)] (2) 


where W is unity in the observed region and zero else¬ 
where. 

Expanding T as in Eq. (|^) produces a set of correlated 
Gaussian variates dim for the signal and un im for the 
noise. These combine into power spectrum coefficients 


Cl 


1 

21 + 1 
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whose statistical properties differ from the ones of the C;. 
We therefore refer to these quantities as pseudo-Ci. In 
what follows we will discuss the statistical properties of 
these quantities, restricting ourselves to a presentation of 
results. We relegate detailed derivations and implemen- 
tational issues to a future publication M . 

If we assume white noise, = C/Gthen each term 
in the sum Eq. ( 0 ) has expectation value 


2 _ lll'm' Cl' \Wl,m' lm\^ + C^f^ djW^ 

21 + 1 


P) 


(4) 


where we abbreviate the polar part of the Vim as Aim = 
\J ^ Wi'm' Im is the matrix element of IT 

in Eq. (|^) in a spherical harmonic basis and integrates 
over the observed region of the sky. Note that the cross¬ 
term which comes from expanding the square in Eq. (|^) 
has vanishing expectation value because signal and noise 
are assumed to be uncorrelated. 

To elucidate the correlation structure between the 
dim we omit the gn im for simplicity and write Eq. (H) 
as a quadratic form in the independent normal vari¬ 
ates aim = GimlVC- This defines the coupling ma¬ 
trix M-d) such that Ci = a*Md)a. Eor the full sky 

l2m2 ~ (2I+1J SmIn tho CaSe 

of azimuthal symmetry is block diagonal: 

(5) 

where = Wi'm imVCv- Each block is the Carte¬ 
sian product of V with itself. Moreover, if Wn has the 
same azimuthal symmetry as IT, this argument can be 
extended to include the noise. Eq. ||is the key fact which 
allows the derivation and cheap evaluation of the exact 
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results we obtain. (To compute the formulas we present 
in this Letter, generating the IT matrices is the most 
costly operation. For a maximum I of 1024 this takes 
just over I minute on a single RIOOOO CPU.) 

It follows that while the dim + djsr im terms in Eq. ^ 
are correlated for different I, they will be uncorrelated 
for different m. As a consequence, we can view the Ci as 
sums of independent Chi-squared (y^) variates, each with 
one degree of freedom but different expectations. This 
allows us to solve for their statistical properties exactly. 

We can compute the distribution of these sums an¬ 
alytically using the method of characteristic functions. 
Defining Sm = (2lcr^m)y produce a closed form solu¬ 
tion in terms of incomplete gamma functions 7 ( 0 , 0 ;): 


P(Ci) = A(') ^ 




(i,(s«-s«)a) 




( 6 ) 


where the primed product symbol only multiplies fac¬ 
tors which have m ^ m' and the normalization constant 

i, ^(,) = Tf nL, ■'£ 

IS ^ — 2 r(f) 

We can then obtain the cumulants of the pseudo-C; as 
Kn = 2 "“^(n — 1 )! moment can be writ¬ 

ten in terms of these cumulants. We give the following ex¬ 
pressions for the mean, variance, skewness 81 = I 

<(AC,)U^ 

and kurtosis (32 = as examples: 

((ZaCz) ) 


.4 

m ^ Im 


{Ci) = Em^L^ ((ag)2) = 2E^^; 


/ 3 i =22 




/32 = 12^ 


(7) 


The covariance between the pseudo-C; can be written 
as (^ACiACr^ = trMd'>MC'> which reduces to 


{ACiACr ) = E E (8) 

m \ li ) 

in the azimuthally symmetric case. All higher order joint 
moments are computed similarly in terms of traces of 
products of M. of various I ||^ . 

There are some important situations where the noise 
pattern does not follow the azimuthal symmetry of the 
survey geometry. In the case of the Planck satellite the 
scanning strategy is approximately centered on the eclip¬ 
tic poles, while the Galactic cut is tilted through « 60° 
with respect to this. In this case Eq. (^) becomes an ap¬ 
proximation. We found it to be very accurate indeed to 
continue using these distributions with the af^ computed 
for an asymmetric ITv, even for a strongly asymmetric 
noise pattern. This approximation will be worst in the 
least interesting, noise dominated regime at very high 1 . 
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Note that the and hence the {Ci) remain exact (be¬ 
cause the Aim are independent of the azimuth) but the 
remarks leading to Eq. (^) are no longer exactly true. For 
applications the final justification comes from the excel¬ 
lent agreement we find when we check against our Monte 
Carlo simulations. 

To test our results we performed 3328 Monte Carlo 
(MC) simulations for a high resolution CMB satellite, 
such as MAP or Planck (resulting in a sky coverage com¬ 
parable to COBE). We simulated realizations of the CMB 
sky in the standard cold dark matter model(nm = 1, 
Vlbh? = 0.015, Hq = 70km/s/Mpc). From these maps 
we carved out a ±20° Galactic cut and contaminated the 
remaining area with spatially modulated Gaussian white 
noise of maximum rms temperature 124^A' per pixel of 
characteristic size 3.4 arcminutes. We use a tilted noise 
template Wn = sin 9 e, where 0 b is the ecliptic lati¬ 
tude, as a simple model of the noise pattern which would 
result from scanning along meridians through the ecliptic 
poles. We then Fourier analyzed these maps and stored 
the resulting Ci. A Kolmogorov-Smirnov test failed to 
detect deviations between the distributions of this MC 
population and Eq. (^ at 99% conhdence, which vali¬ 
dates our semi-analytical expressions. 

To give a visual impression of the resulting probability 
densities we show four of them in Figure ^together with 
the results from the MC simulations. Also shown in this 
Figure are the distributions which the Ci would follow 
in the full sky case as well as the commonly used Gaus¬ 
sian approximation. These are mean adjusted to account 
for the lost solid angle due to the Galactic cut. At Z < 30 
the difference is striking. For higher I the Gaussian ap¬ 
proximation becomes better as higher moments die away 
by dint of the Central Limit Theorem, but there remain 
visible systematic differences to the true distributions. In 
particular, there is a residual shift in the mean and the 
approximations tend to be slightly narrower than the his¬ 
tograms for very high I « 1000. 

This becomes a more quantitative observation when 
looking at the percentage discrepancies between the 
mean and variances as a function of I in Figure |^. The 
discrepancy is of the order of 1 % in the mean and 5% in 
the standard deviation on most scales, except for I < 70 
where the effect is larger. These discrepancies are im¬ 
portant at the level of precision of future almost full sky 
missions. For medium and small sky coverage the mode 
couplings are stronger and we expect this to have an even 
larger effect on the probability distributions. We also 
compare the skewness /3i and kurtosis P 2 of the dis¬ 
tributions to our distributions. The percentage difference 
is larger than for the first two moments but arguably less 
important at large I, since /3i and (32 decay as {21 + 1 )“ 2 
and (2Z ± 1)“^, respectively. 

As a first application we study the effect of approx¬ 
imating the likelihood for parameter estimation. Since 
our distributions have the correct means, we simply mul- 
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FIG. 1. The Pseudo-Gi distributions, Eq. ^ for the stan¬ 
dard cold dark matter model (shaded), the (solid lines) and 
the Gaussian (dashed) approximations compared to Monte 
Carlo simulations (histograms) for 1=2,200,310,1000. 


tiply them together for a simple, unbiased approximation 
to the likelihood. This is conservative since the marginal 
distributions have all correlations integrated out and we 
will therefore overestimate the error bars on the G/. Us¬ 
ing this likelihood, as well as the Gaussian and approx¬ 
imations, we attempt to estimate the baryon parameter 
Uf, (holding all other parameters constant) from several 
randomly selected realizations in our MG pool. The re¬ 
sults are shown in Figure |[ As expected, Gauss and 
consistently hnd estimates which are biased about 1.6% 
high, 3 standard errors of the mean away from the true 
value, while our likelihood gives a perfect fit. Since the 
moment discrepancies depend on the underlying cosmo¬ 
logical theory, this level of bias is only indicative of the 
general level of error introduced by using the Gaussian 
or approximations. 

We note that numerical techniques have recently been 
developed |11| which allow the computational solution 
of the power spectrum estimation for high I applications. 
These methods are efficient for the case of almost full sky 
observations. A purely numerical approach still requires 
signihcant computational resources, especially if there is 
signal in modes with I > 1000. The results presented 
in this Letter should be seen as complementary to such 
calculations. An analytical framework admits a more 
fundamental approach to understanding and is a useful 
yardstick against which numerical work can be tested or 
from which approximate methods can be derived. 

To illustrate, we suggest the following computationally 
cheap and accurate approximation recipe, motivated by 
the fact that higher moments of the Ci distributions die 
away quickly with higher 1 : on measuring the Ci from 
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FIG. 2. The solid lines are, top to bottom, the percent¬ 
age discrepancy between the mean, standard deviation, skew¬ 
ness, and kurtosis of the approximation and the Pseudo-C* 
distributions. The stars are the {Ci) computed from 3328 
Monte Carlo runs, showing excellent agreement only limited 
by Monte Carlo noise to better than 0.1 % for all Z. 


FIG. 3. One parameter likelihood estimates of fit- Shown 
are the means (solid) within ±3 standard errors in the mean 
(dashed) of estimations from 17 realizations (diamonds) of the 
sky. Both the x^ (left) and the Gaussian (right) approxima¬ 
tions produce a bias, while our approximation (center) has no 
detectable bias. 


the sky, fit a smooth curve through them and use this 
fit as the input theory for our framework. Then calcu¬ 
late the corrections to the sample means and variances 
using Eqs. (||j^). These can in turn be used for simple 
X^ fitting in the usual way, at least for high I > 100. 
Since the discrepancies are small for large sky coverage 
this will produce a fit which is accurate to second order 
in the discrepancy, safely within the regime of accuracy 
envisaged for modern satellite missions. 

Apart from the obvious applications to experiments 
with small and medium sky coverage such as balloons or 
ground-based missions, many further uses of this frame¬ 
work are conceivable. For example, one could 1) extend 
this treatment to multi-parameter fits, 2) design ’opti¬ 
mal’ scanning strategies, encoded in Wn, 3) use Eq. 
once a theory is estimated to check I by I for consistency 
with the assumption of Gaussian primordial fluctuations 
and 4) assess more realistically if secondary anisotropies 
will be detectable with future CMB missions. Finally, 
all ingredients are there to refine the approximation to 
the joint likelihood we used in this Letter by taking into 
account the covariances Eq. (||), for example by using 
a multivariate Edgeworth expansion around the peak of 
the likelihood. 

To summarize, we have presented a unified theoret¬ 
ical framework for the study of power spectrum de¬ 
terminations of balloon, ground based and satellite ex¬ 
periments. We go beyond current approximations and 
present a semi-analytic formalism for the computation of 
sampling distributions of the Ci for any Gaussian cos¬ 
mological model and a large and important class of sur¬ 
veying strategies. We show that applying this method 
to the estimation of fZf, from simulated data is unbiased 
and hence superior to commonly used analytic approx¬ 


imations which bias the result at the percent level. A 
number of applications and extensions of this formalism 
remain to be explored |]l^ . 
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